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We consider stochastic processes where randomly chosen particles with positive quantities x, y (> 
0) interact and exchange the quantities asymmetrically by the rule x' — c{(l — a)x + by}, y' = 
d{ax + (1 — b)y} (x > y), where (0 <) a, b (< 1) and c, d (> 0) are interaction parameters. 
Non-integer power-law tails in the probability distribution function(PDF) of scaled quantities are 
analyzed in a similar way as in inelastic Maxwell models(IMM). A transcendental equation to 
determine the growth rate 7 of the processes and the exponent s of the tails is derived formally 
from moment equations in Fourier space. In the case c = d or a + b — 1 (a / 0,1), the first-order 
moment equation admits a closed form solution and 7 and s are calculated analytically from the 
transcendental equation. It becomes evident that at c = d, exchange rate b of small quantities 
is irrelevant to power- law tails. In the case c 7^ d and a + b ^ 1, a closed form solution of the 
first-order moment equation cannot be obtained because of asymmetry of interactions. However, 
the moment equation for a singular term formally forms a closed solution and possibility for the 
presence of power-law tails is shown. Continuity of the exponent s with respect to parameters 
a, b, c, d is discussed. Then numerical simulations are carried out and campared with the theory. 
Good agreement is achieved for both 7 and s. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.20.Dd, 89.65.-s 



I. INTRODUCTION 

A mechanism of power-law emergence has been taken 
a strong interest for many researchers because of its ubiq- 
uity in nature. A number of attempts have been devel- 
oped to explain the advent of various power-law distri- 
butions P, Qi 0, 3 , but underlying physics has not been 
clarified yet. 

Recently, inelastic Maxwell models (IMM) [HS0 have 
been studied extensively and provides a new mecha- 
nism of power laws. In IMM, randomly chosen parti- 
cles undergo binary inelastic collisions and non-integer 
power-law tails appear in probability distribution func- 
tion(PDF) of particle velocities. This processes are de- 
scribed by the Boltzmann equation with a velocity inde- 
pendent collision rate and are analytically tractable. The 
moment equations of IMM form a closed set which is able 
to be solved sequentially as an initial value problem. A 
nontrivial power-law exponent which is the function of a 
dissipation parameter is determined from a transcenden- 
tal equation. This power law differs from usual critical 
phenomena in one major point that the system has a 
power-law region rather than a ciritical point. Therefore 
fine-tuning of parameters is not necessary and the power 
law is easy to be observed inside the region everywhere. 
ben-Avraham et al. extended IMM to a symmetrically 
linear collision rule. 

On the other hand, power laws are widely observed in 
social and economic phenomena fToL ITU hj. It is well- 
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known that wealth distributions such as capitals and in- 
comes obey a power law in high-wealth range, which is 
called Pareto's law® H Q tH 13. It is also recog- 
nized that the size 5 of cities satisfies power law 1/5 
in the cumulative distribution function, that is, Zipf's 
lawfTol IttI ITc|. A multiplicatively interactin g st ochastic 
process is one candidate to explain the laws [17], [n|- In 
economic phenomena, the Matthew effect (the rich gets 
richer and the poor gets poorer) plays an important role, 
which results in asymmetry in interactions. Ispolatov et 
al. |20| have discovered that a power-law distribution 
with the exponent of unity in PDF arises in multiplica- 
tive processes of greedy exchange, where the rich always 
gets richer. The purpose of this paper is to study a model 
of asymmetrically multiplicative interactions and to clar- 
ify the influence of asymmetry. Note that the model in- 
cludes IMM 0, H 0> symmetirically multiplicative in- 
teraction(SMI) model [a , and greedy multiplicative ex- 
change(GME) model [2(j as special cases. 

The paper is organized as follows. In Sec[nj we begin 
with introducing the model. The Fourier transform of 
the master equation is performed. PDF is assumed to be 
sum of the regular term and the singular term. Then, 
a transcendental equation is derived from the singular 
terms. In Sec. IIIII the growth rate 7 of the processes and 
the power-law exponent s of tails are discussed in three 
cases I, F, and II. In case I and F, 7 and s are com- 
puted explicitly by solving the transcendental equation. 
In case II, moment equations do not form a closed set, 
and 7 and s cannot be calculated analytically. Thus nu- 
merical simulations are performed. The good agreement 
between the theory and simulations is achieved. In the 
last section, we discuss our results. 
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z - cbx - c(l-a)y — 



z - d(l-b)x - da y = 




z - c(l-a)x - cby — 



FIG. 1: Integral lines in Eq. @. 



II. TRANSCENDENTAL EQUATION 

Let us consider the stochastic processes that distin- 
guishes two particles in the manner of the magnitude 
of quantity asymmetrically: When a particle of posi- 
tive quantity x (> 0) interact with a particle of quantity 
y (> 0), post-interaction quantities x' and y' are given 
by 



x' \ _ f c(l - a) cb 
y' ) ~ I da d(l - b) 



(*>V) (i) 



where < a, b < 1 and c, d > are interaction parame- 
ters representing amplification rates and exchange rates 
of larger and smaller quantities, respectively. Two parti- 
cles are selected randomly. There are two trivial cases (i) 
a = b = or a = b = 1, (ii) a = 0, b = 1 or a = 1, b = 0. 
In the former case, two particles experience no interac- 
tions and PDF becomes log-normal. In the latter, one 
particle ends up taking all the quantity and the others 
having null quantities. Hereafter we omit these cases. 

The normalized distribution function f(z,t) obeys the 
master equation. 



df(z,t) 
dt 



+ f(z,t) = 



dy / dxf(x,t)f(y,t) 
o J y 

x [ 6(z — (c(l — a)x + cby)) 
+5(z - (dax + d{l - b)y)) ] .(2) 



The integration is performed along the solid line illus- 
trated in Fig. ^ The kink at y = x stems from asym- 
metry of the model. The Fourier transform g(k, t) — 
J °° dze tkz f(z,t) ofEq. © is performed by changing vari- 
ables p = x — y, q = y and using the formula 



dpe lkp = 7T(5(fc) + i v.p. 



where v.p. means the Cauchy principal value. It follows 
that 



—g(k,t)+g(k,t) 



-[g(c(l - a)k,t)g(cbk,t) + g(dak,t)g(d(l - b)k,t)} 



v.p. 



2n 



v.p. 



dk' 
dk' 



1 



k' - c(l - a)k 
1 



g{k', t)g(c(l -a + b)k- k', t) 



k 1 — dak 



g(k',t)g(d(a + l-b)k- k',t). 



(3) 



The v.p. terms in r.h.s. of Eq come from the inte- where m n (t) = dzz n f(z,t) is the n-th order moment, 
gration of kinked lines, and represent asymmetry of the Even though the moment equations do not form a closed 



model. Next, g(k,t) is expanded as 



set generally, they can be derived formally as 



g(k,t) = l + 



(ik) n 



(4) 



dt 



m n (t) - X n (t) m n (t) 



1 71—1 

1=1 



{c"(l 



d n a l (l - b) n - l }mi(t)m n ^(t) (n > 1), 



(5) 



*»(*) = ^[c n {(l^a) n + b n } + d n {a n + (l-b) n }}-l 
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n ^) E ( " ) W - a ) V ~' + ^ - V -P- J" dk 'y9 {l) (k', t)g {n - l \-k', t) (6) 



where g^ l \k, t) denotes the I— th order deirivative In real x-y space, the pseudo eigenvalue X n (t) is expressed 



as 



ffW(fe,*) = < , E^r m '+«(*)- ( ? ) 

' n! 

n=0 

I 



K(t) = - [c n {(1 - a) n + b n ) + d n {a n + (1 - &)"}] - 1 



1 1 

2 TO 



l n n—l n I svJi—l 



y l x n - l )f(x,t)f(y,t). (8) 



Since f(x,t) includes information of all-order of mo- 
ments, the moment equations do not form a closed set in 
the presence of the integral terms of Eq. JSJ . In IMM and 
SMI, however, the integral terms vanish. Therefore, A„ 
is constant and the moment equations (JSJ) form a closed 
set. 
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We are interested in similarity solutions of the form 



1 s(M) 



$(r?) 



(9) 



where 7 is the growth rate (scaling parameter) of the 
system, and £ = ze~ 7 * and 77 = ke 7t are scaled variables. 
The scaled PDF $(77) satisfies 



+ *(»7) 



[$(c(l - a)rj)${cbri) + &(da<n)$(d(l - 6)77)] 



2tt 
'2tt 



dfc' 

1 

dk' 



1 



77' — c(l — a)rj 
1 



$(r/)<I>(c(l - a + b)r] - ?/) 



77' — darj 



$(?/)$(d(a + 1 -6)77-77')- 



(10) 



Here we assume that the function $(77) is described by 
the sum of the regular and singular components 

$(77) = ^regulariv) + $ singular (f?) , (H) 
^2, («77)™ 

$regniar(f?) = 1 + V- L - T -^„, (12) 

n=l 

(7TS \ 
-i — sgn(rj)j T(-s)\i]\ s (13) 



where /_t n is the 71-th order moment of $(77), s a 
non-integer exponent, C the normalization constant, 
r(— s) (s > 0) the gamma function, and sgn(rj) the sig- 
nature of 77. The form of the singular component Eq. i|13|) 
is given by the Fourier transform of = C/£ 1+s . 

The leading small-77 behavior of the singular component 
& singular (v) ~ M s reflects the tail of the scaled PDF 
~ l/£ 1+s as £ -> 00. Substituting Eq. JTTJl into 
Eq. I|10(l . we obtain relations: 
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Ofa)-> 7 = Ai= i[c(l-o + 6) + d(o + l-6)] + i(c-d)(l-o-6)A-l, 



7 s = A s = {c(l - a)} s + {da} s - 1 



where 



/>oo />OC 

/ d& / da(Ci-6)*(a)*(6 

Jo J{ 2 



A= — 
Mi Jo 

and the eigenvalue A„ is given by 

A„ = \ [c n {(1 - a) n + b n } + d n {a n + (1-6)"}]-! 



(0< A<1), 



1 1 

2 /U n 



(14) 



n-l , s 

0(r? n ) -> (rvy ~ A n )/U„ = - ]T " {c"(l - a)<6"- ; + dV(l - 6)""'} Mi /z„_; (n > 2), (15) 

1=1 ^ ' 



(16) 



(17) 



In the process of deriving Eq. I|ltj|) . we use a convergence condition; $(r/) — > — > oo) and the formulae to calculate 
regular-singular integrals 



u.p. / drf'—j— |j/| s $ resu ;ar((a + 0)v - v') = ta.n( — ) tt sgn(arj) \m]\ s $ regu i ar (/3r]), 

J-oo V a7 l * 
° 1 7TS 

u.p. / dri'—j— sgn(ri') \rj'\ s $ regu i ar ((a + ft)?] - 7]') = - cot(-) n \ar]\ s $ regu i ar (f3r]). 

> ?7 ar l £ 
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In the caluculation of deriving moment equations, 
singular-singular integrals are able to be neglected be- 
cause the higer-order singular terms in Fourier space are 
less singular in real space. That is, it is possible to lin- 
earize the moment equations for the singular term, which 
is the key ingredient of derivation. Notice that moment 
equations l|14fl and l|15fl for regular parts can be directly 
derived from the master equation J5J without recourse to 
Fourier transform. 

The transcendental equation of the model with asym- 
metrically multiplicative interactions is formally ob- 
tained from Eqs. (|14|l and (|16|) by eliminating 7 



\ x s = {c(l - a)} s + {da} s 



s > 1) 



(20) 



When Eq. I|2(J|) has a non-trivial solution in < s < 1, 
the singular component dominates over the regular one 
because of the divergence of all the moments. Therefore 
the parameter 7 is not determined by Eq. (|14f> . As origi- 
nally reported by ben-Avraham et instead, 7 takes 
the minimum value, which is realized when the line 7s 
comes into contact with the curve {c(l — a)} s + {da} s — 1. 
Consequently, the transcendental equation in the case 
< s < 1 becomes 



{c(l -a)} s ln 



Ml -a)}' 



{da} s In 



{da}' 
(0 < s < 1). (21) 



III. 



GROWTH RATE AND EXPONENT 



In this section, the growth rate 7 of the processes and 
the exponent s of power-law tails are investigated. Three 
cases I, I', and II are discussed separately in accordance 
with the type of the transcendental equation. Equa- 
tion 114f) tells us that the first-order moment equation 
admits a closed form solution and the first-order eigen- 
value Ai becomes constant in the case c = d (I) and 
a + b = 1 (o ^ 0, 1) (I'), while it does not otherwise (II). 
When Ai is constant, 7 and s can be calculated analyt- 
ically, which are done in subsections A and B. Case II 
is treated in subsection C. In this case, explicit calcu- 
lation of 7 and s becomes impossible at s > 1. When 
< s < 1, on the contrary, the selection of the minimum 
growth rate leads to the irrelevance of Ai in determining 
7 and s. In this situation, therefore, the transcendental 
equation is given by Eq. i|21|) and 7 and s are computed 
analytically. First, the continuity of s at s = 1 is dis- 
cussed. Then, numerical simulations are carried out and 
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FIG. 2: Phase diagram of the exponent s in the case c — d. 
Solid curves denote contours of s. Power- law tails disappear 
in the hatched region surrounded by dashed curves. 



compared with the theory. 



zero |20(. In the limit a — > or a — > 1, s also goes to zero 
at c ^ 1, which suggests that the cases a = and a = 1 
for arbitrary values of c belong to the same universality 
class as GME. In Appendix, we show this analytically. 



B. Case I' : a + b = 1 (o ^ 0, 1) 

At a + b = 1, Ai = c(l — a) + <ia — 1 and the transcen- 
dental equation becomes 



{c(l - a) + da - l}s = {c(l - a)} s + {tfa} s - 1 (24) 

for s > 1. The equation for < s < 1 is given by 
Eq. (|21[) . Cross-sectional phase diagrams at fixed d for 
s are plotted in Fig. [3J Boundaries between the regions 
with and without power-law tails are described by three 
curves 



1 



1 



1 — da 
1-a ' 



(25) 



Therefore the area without power-law tails is unbounded 
when d < 1, while it is bounded when d > 1. Asym- 
metry between c and d gives rise to asymmetry of phase 
diagrams. 



Case I : c = d 



Case II : c d and a + b 7^ 1 



When c = d, first-order eigenvalue Ai is given by c — 1. 
Then, the transcendental equation reads 

(c - l)s = {c(l - a)} s + {m} s - 1 (s > 1), (22) 
{c(l -a)} s ln 



-j— — + {ca} s In 

|c(l — a)}"/ \{ca} 5 

= 1 (0 < s < 1). (23) 



It should be emphasized that the transcendental equa- 
tion, the growth rate 7, and the exponent s are indepen- 
dent of b, exchange rate of smaller quantities. When am- 
plification rates of larger and smaller quantities arc same 
(c = d), only exchange of larger quantities determine 
similarity properties of processes. Although Eqs. I|22|) 
and l|23|) are formally same as those for symmetric in- 
teractions (c = d, a = b) |8|] , the physical meaning of the 
parameters is completely different. A phase diagram for 
s is illustrated in Fig. [3 Power-law tails exist outside the 
hatched region. Inside the region, s diverges and power- 
law tails disappear. The exponent s varies continuously 
as a function of parameters a, c. The diagram is sym- 
metric with respect to the line a — 0.5. The two points 
c = 1 , a = and c = a = 1 are singular points where 
the transcendental equations become identities. Actu- 
ally, these points correspond to greedy multiplicative ex- 
change(GME) processes and the exponent s is given by 



In this case, the first-order moment equation 114f) does 
not allow a closed form solution. However a transcen- 
dental equation can be derived formally and possibility 
for the presence of power-law tails can be shown. The 
results are given by Eqs. (|14|) and (JTSJ) . The value of 
A in Eq. I|17() cannot be obtained analytically. Therfore 
7 and s are unable to be calculated at s > 1. When 
< s < 1 and the selection of the minimum growth rate 
takes place, all of Aj (i > 1) are irrelevant and 7 and s 
are determined from Eq. I|21|l . 

First we discuss the continuity of s(a, b,c,d) at s = 1. 
Suppose that the power-law tail dominates in determin- 
ing A and we set <F(£) = (1 + e)L 1+ 7£ 2+£ (s = 1 + e), 
where is a normalized power-law PDF and L is a 

lower cutoff in order to avoid the divergence of the inte- 
gral. We obtain 



A S 



1 



l + 2e 



1 (e-0). 



(26) 



Thus 



c(l 



da — 1 . 



(27) 
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This fact indicates s is continuous and satisfies Eq. i|16|) 
identically at s = 1. 

In the general case, < A < 1 and the transcendental 
equation has two nontrivial solutions, Si < 1 and 
S2 > 1 generally. Ernst and Brito [6j considered the 
case and argued that each solution has a different role. 
That is, si is a expansion variable and S2 controls the 
singularity of the tail, 



-fl n + (J, S2 7] S 



(28) 



We find one more reason to explain why S2 determines 
the power-law tail. In the limit c — > d ot a + b — ► 1, 
si — ► 1 and S2 — > s. Assumption that s is a continuous 
function of parameters a, 6, c, d leads to the conclusion 
that S2 corresponds to a nontrivial solution in the cases 
I and I'. 

Now we compare the theory and numerical simula- 
tions. The value of A is computed numerically. Then 
the growth rate 7 and the exponent S2 (and si) are cal- 
culated via Eqs. (|14|> and (|16|) . Obtained results are com- 
pared with those of simulations in Figs. l4l5l We find that 
Eqs. JUJ) and JEJ hold well. 
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FIG. 4: Growth rates 7 calculated through Eq. (11 H and those 
by simulations at c = 1.5 and d = 0.5. Number of particles 
in simulations is N = 10 6 and number of interactions is T = 
50iV. 



IV. SUMMARY AND DISCUSSIONS 

In this work, we have investigated asymmetrically mul- 
tiplicative interactions (AMI) processes analytically and 
numerically. The model includes greedy multiplicative 
exchange (GME) model |20j, inelastic Maxwell model 



(IMM) HUB], an d symmetrically multiplicative inter- 
actions (SMI) processes [g| as special cases. The relations 
are summarized in Table [I] 

In the processes of deriving the transcendental equa- 
tion, we encounter the fact that the first-order moment 



7 




logio® logio© 

FIG. 5: Double logarithmic plot of versus £ at a — 0.05 (a), and a = 0.2 (b). Other parameters are given by b = 0.9, 

c = 1.5, and d = 0.5. Number of particles is TV = 10 6 and number of interactions is T — 50N. Dashed and dash-dotted lines 
are those of the slope S2 and sj, which are determined from the transcendental equation. 



TABLE I: Relation between AMI and other models. 



Model 


Interaction parameters 


GME 


c = d = 1, a — 


IMM 


c — d = 1, a = b 


SMI 


c = d, a — b 



equation does not admit a closed form solution generally. 
However, it is possible to linearize singular terms. Thus 
we can formally construct the transcendental equation 
and show the possibility for the existence of non-integer 
power-law tails in PDF. It becomes evident that when 
c = d, the exchange parameter b of small quantities is 
irrelevant to power-law tails. In the general case c =/= d 
and a + b =/= 1, the growth rate 7, the exponent s, and 
PDF are determined consistently via Eqs. (|14fl and (|16|l . 
The exponent s turns out to be a continuous function of 
parameters a, 6, c, d. 

This model is applicable to many kinds of asymmetri- 
cally interacting processes such as particle systems and 
biological phenomena. The most important would be 
wealth distributions in economic systems, where the 
Matthew effect gives rise to asymmetry of interactions. 
The model is quite general if terms are replaced as " parti- 
cle" — ► " agent" , " quantity" — * " capital" or " income" , and 
"interact" — > "trade" or "deal". Power laws are widely 
observed in economic phenomena 1 1 ll . Some of them 
might be explained by this model. 



APPENDIX: GENERAL GREEDY EXCHANGE 
PROCESSES 

Ispolatov et aZ.[2(J have reported that greedy multi- 
plicative exchange processes where a = and the rich 
always gets richer exhibit a power-law behavior 1/z in 
PDF. They refered only to the case c = d = 1, but as 
discussed in the case I, the tail 1/z with s = is also re- 
alized at c 7^ 1 and d 1 . Here we extend their approach 
to general greedy exchange processes a = 0, 1, d ^ 1 
and examine the condition for f(z,t) oc 1/z. The master 
equation is 



df(z,t) 
dt 



= ~f(z,t) 



(A.l) 



1_ [— 

cb 



1 



d(l - b) 



[ Z C{ ] b a)X ..t 



dxf(x,t)f I ^ aX , ,t 



d(a+l-&) 

In the limit of a — > 0, Eq. (|A.1(I reduces to 



d(l -by 



df(z,t) 
dt 



-f(z,t) 

'7b (r^T' 1 



(A.2) 



+ 



1 



d(l-b) J \d(l-b) 



-,t 



dxf(x, t). 



d(i-b) 
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Performing the transformation (z — cx)/cb — > x and using 
the normalization condition J Q dxf(x,t) = 1, we have 



df(z,t) 

at 



dxf(x, t) 



7M) 



dxf(x, t) 



+ 



-f 



d{l-b) J \d(l-b) 



(A.3) 
dxf(x, t). 



In the limit a — > 0, divergent terms vanishes and the 
solution of Eq. HA.3|) is obtained by 



D 



D 



-1 



In 



rf(l-b) 



(A.4) 



Equation l|A.4|) shows that the condition d(l — b)/c < 1 is 
necessary to ensure f(z, t) > . It becomes evident that 
the condition that general greedy exchange processes has 
the distribution 1/z is given by 



d{\ -b) < c. 



(A.5) 



We find that the condition l|A.5[) is always satisfied at 
c = d. 



Next, we estimate the lower and upper cutoff's z±, Z2 of 
the distribution Eq. (|A.4|I in order to maintain the nor- 
malization condition of PDF. The moments are evaluated 



as 



m (t) ~£ a <fe/M)~ J ln (f 1 1 



mil 



~ / z ? dzzf(z,t) ~ y(z 2 - Zl) = mi(0)e 
Using Z2 >> zi, we get 



7* 



-21 



mi(0) H-i/DH 
D 



*~^te« (A.6) 



We have carried out numerical simulations of general 
greedy exchange processes. The results are shown in 
Fig. EJ The distribuitons oc 1/^ is observed for 

several tens order of magnitude. It should be empha- 
sized that in the case c > 1, the tails with s > appear 
asymptotically (£ — > oo). With increasing c, s decreases 
and fluctuations grow. This indicates that the tails come 
from fluctuation due to amplification of large quantities. 
However, detailed analysis of the tails remains for a fu- 
ture work. 
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FIG. 6: Double logarithmic plot of versus £ at a — 0.0, 

b = 0.9, and d = 0.5. Number of particles is TV = 10 6 and 
number of interactions is T = 50N. 



